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We analyze wave propagation in coupled planar waveguides, pointing specific attention to modal 
cross-talk and out-of-plane scattering in quasi-planar photonics. An algorithm capable of accurate 
numerical computation of wave coupling in arrays of planar structures is developed and illustrated 
on several examples of plasmonic and volumetric waveguides. An analytical approach to reduce or 
completely eliminate scattering and modal cross-talk in planar waveguides with anisotropic materials 
is also presented. 



I. INTRODUCTION 

On-chip communications, surface plasmon optics, and 
Si photonics are all examples of planar optics, where the 
optical radiation is controllably guided on the plane of 
a photonic chip. A number of planar optical elements 
including lenses, mirrors, and on-chip waveguides, both 
plasmonic and dielectric have been recently designed, 
fabricated, and characterized^ ^^^i 10 ' 11 ' 1 ^ . How- 
ever, with a few exception o 11 ' 13 ' 14 , the majority of recent 
studies focuses on in-plane propagation of light and ne- 
glects out-of-plane scattering of radiation. In this work 
we analyze the out-of-plane light scattering and modal 
cross-talk due to effective index change inside planar 
waveguides and demonstrate that this scattering may 
substantially affect the propagation of confined modes 
in complex planar systems. We design a numerical ap- 
proach to solve the problem of scattering and modal 
cross-talk in planar or quasi-planar structures that con- 
tain anisotropic elements, and present the technique to 
reduce or completely eliminate the scattering and cross- 
talk with anisotropic metamaterials 1 ^. 

The problem of out-of-plane scattering has been of con- 
sistent interest to photonic community. Although a num- 
ber of finite-difference and finite-elements techniques, 
available today, can successfully solve the problem of 
scattering in relatively small geometries (< lOAo, with 
Ao being free-space wavelength) , analysis of wave propa- 
gation in an extended system is beyond the capabilities 
of methods that rely on finite-size meshing of space/time. 
One of the ways to reduce memory requirements to cal- 
culate the field in an extended structure with moderate 
number of scattering interfaces is to implement some sort 
of wave-matching technique where the modal spectrum 
is constructed to satisfy the solutions of Maxwell's equa- 
tions in the space, and only boundary conditions at scat- 
tering interfaces are enforced, resulting in calculations of 
amplitudes of the modes. Effectively, modal expansion 
can replace the need to calculate all field components at 
every point of space with the need to calculate modal 
amplitudes in every region of space. 

One of the first descriptions of wave-matching ap- 
proach and its applications for highly- conductive plas- 
monic guides can be found in Refs i 13 i 16 . Ref3 also de- 
scribes scattering by planar guides with highly symmetric 



crossections. An approach to calculate to modal cross- 
talk and scattering in ID guides was developed in Refs; 17 . 
Scattering by periodically corrugated systems has been 
analyzed with rigorous coupled wave analysis (RCWA) 
inl£. Recently, the generalization of field expansion to 
calculate scattering in plasmonic planar guides has been 
presented inM. However, while the mode-matching cal- 
culations were proven to be highly efficient, the technique 
had failed in the proximity to plasmon resonance condi- 
tion, when the field of a surface wave is highly confined 
to the proximity of metal interface. Here we present a 
wave-matching technique that is capable of solving for 
wave scattering in complex systems formed by coupled 
planar waveguides. 

The rest of the manuscript is organized as follows. In 
Section 2, we present the mode structure of an arbitrary 
planar guide used in this work (FigHJ. Section 3 is de- 
voted to the development of mode-matching technique in 
quasi-planar system comprising a uniform in y-direction 
array of planar guides. The presented numerical ap- 
proach is illustrated on examples of light propagation in 
several plasmonic and metamaterial systems in Section 
4. Finally, Section 5 develops the formalism of truly pla- 
nar photonics where out-of-plane scattering and modal 
cross-talk are not possible and presents an approach to 
design extremely-low-scattering plasmonic circuits. 

In this work we use the following notations for the elec- 
tric and magnetic fields in the system: the total electric 
(E) and magnetic (H) fields are shown in italic font; the 
fields of modal components contributing to the total field 
are represented with symbols E and H; the components 
of an individual mode in a particular layer of the mul- 
tilayer stack are represented with calligraphic symbols £ 
and TL. 



II. MODAL SPECTRUM OF PLANAR GUIDES 

We start from analysis of modal spectrum of a planar 
waveguide, schematically shown in Figflja). The struc- 
ture comprises a set of N planar layers with layer inter- 
faces parallel to yz plane, with j th layer occupying the 
space between Xj-\ < x < Xj, and having (uniaxial) di- 
electric permittivity described by a diagonal tensor with 
diagonal components ij = {e* x , e^ z , ej z }- We assume 
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FIG. 1: (color online) Schematic geometry of the multi-layered structures and electromagnetic mode types used in the 
manuscript. Geometry of the single multilayer stack is shown in (a); panel (b) explains the composition of top, bottom, 
and guided modes; field profiles in the outside layers of the multilayer are shown; the interface between two multilayer stacks 
is shown in (c) 



that xq = — oo, and xn — oo. In this section we consider 
the structure that is infinitely extended in the yz plane. 

The electromagnetic fields in this layered system 
can be represented as a set of transverse-electric (TE) 
and transverse-magnetic (TM) polarized waves (modes). 
Each mode of the multilayer constitutes a solution of 



Maxwell's equations that is finite for — oo < x < ocJ^. In 
homogeneous layered structure the mode can be parame- 
terized by a combination of (i) its polarization (TE/TM), 
(ii) the in-plane components of the wavevector (k y ,k z ), 
and (iii) a set of layer-specific complex coefficients {af} 
playing the role of amplitudes of the mode components: 



&(k y ,k z ) = < 



a^£ i (k X \i (ky 7 kz) > ky, k z ) + £ \ ( k x -\(k y , k z ), k y , k z ), x < x\ 

Q>j £j{kx;j{k y , kz) , ky , k z ) + a- £j( k x -j(k y ,k z ) 7 k y ,k z ), Xj—i < x <C xj 



M{k y ,k ZJ 



, aJj£ N {k x . N (ky, k z ),k y ,k z ) + a N £ N (-k x . N {k y , k z ), k y , k z ), x N _i < x 
ti^Hi {k x -,i {ky , kg"), ky, kz) + cii Wi ( kx;i(k y , k z ) , ky , k z ) , x < x\ 

o^1~Lj {k X ;j(ky , k z ) , ky , k z ) -t- cij 7~tj{ kx-,j{ky,kz),ky,kz), ^j— i < x < Xj 
> a%HN(k x -N(k y , k z ), k v , k z ) + a^HN(-k x; N (k y , k z ), k y , k z ), xjv-i < x 
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where the fields and dispersion relations in each layer are described by: 



TE — polarized waves 



TM — polarized waves 



{o, 



kz 



-ky} 







{ (ky -+- k z ), k x -jk y , k x -jkz\ 



_ ...t..T 



7^ : 



ej x k z , 



e xx k k 
7^ — J 



(3) 



(4) 



Note that for a given mode only two of the amplitudes transfer matrix method 19 : 

are independent of each other. Indeed, the remain- 
ing amplitudes can be calculated using the well-known 



(l-KMtf (l + K 3 )/^ 



,(5) 
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with polarization-dependent parameters <ftj , ctj , and 
Kj given by (j>f = exp[i(fc x;3+1 ± k x]j )xj], aJ E = 
1/2, a™ = ef/{2e?l x ), Kj E = k x . d /k x;j+u K™ = 

For multilayer systems, layer-specific transfer matri- 
ces can be multiplied together, yielding the transfer- 
matrix relating the fields in any two layers of the mul- 
tilayer stack, and thus solving the problem of reflection 
and transmission of a plane wave by the multilayer com- 
posite. The singular solution that corresponds non-zero 
scattered waves on both sides of the multilayer structure 
(aj~,a^r 7^ 0) with zero incident fields (a^ = = 0) 
corresponds to eigen (guided) mode of the stack^. 

For each polarization, full spectrum of modes sup- 
ported by the stack includes three groups of waves [see 
FiglUb)]. The first group contains a discrete spectrum 
of guided modes, exponentially decaying into first and 
last layers (a^ = a^- = 0, a7/ = 1). Here we charac- 
terize these waves by the in-plane components of their 
wavevectors {k y ,k z }. 

The remaining groups of modes contain the continuum 
of waves, known as open- waveguide mode o 13 ' (bulk 
modes). The first of these groups represents the modes 



originated by a plane wave that is incident on the lay- 
ered structure from the top layer (af = 0, ajf = 1), 
while the second group represents the wave incident on 
the structure from the bottom layer (a x = 1, = 0). 
Here, the modes of the first group ("top modes") are 
parameterized by the real-valued k x; N in the top (N th ) 
layer, while the "bottom modes" are parameterized by 
the real- valued k x -i in the bottom (1 st ) layer. 

In the limit of symmetric distribution of permittivity 
e(x) = e(—x), the proposed here spectrum of top and 
bottom modes is equivalent to the proposed earlier- 
combinations of "standing wave" modes with symmet- 
ric and anti-symmetric x-profiles. However, in contrast 
to the latter, the combination of top and bottom modes 
is more easily generalizable to the case of non-symmetric 
(such as plasmonic) planar guides. Note that in major- 
ity of previous studies of plasmonic structures^, bottom 
modes were explicitly omitted. As explained below, this 
omission becomes crucial in the regime of strong surface 
plasmon polariton (SPP) scattering, e.g. in proximity to 
SPP resonance or in plasmonic step geometry (Figl3|). 

Overall, the field inside the guiding structure can be 
written as: 



E 



E 

+ /o° 
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E(k y , ki q) ) + A<-^E(k v ,-k z q) ) 

(k x )E(ky,k z (k x , N )) + A^' (k x )E(ky,-k z (k x , N )) 
+ A ht ^ + (k X )E(ky,k Z (k X , 1 )) +A ht ^(k X )E(ky,-k Z (k X . 1 )) 

A^ + M(k y , k z q) ) + Ato'Wfo, -ki q) ) 
A to P + (fc x )H(fc y ,fc z (fc x;JV ))+v4 to P"(fc x )H(fc y ,-fc z (fc x;JV )) 

+ A htm+ (k x )U(ky,k z (k x . 1 )) +A ht ™~ (k X )W(ky,-k Z {k x;1 ))\ ilk, 



(6) 



Here we assume that all the modes in the layered 
material have the same value of k y . This assumption 
does not limit the applicability of the developed tech- 
nique since due to translational symmetry, any solu- 
tion of Maxwell's equations in the coupled waveguide 
set can be represented as a linear combination of so- 
lutions for specific values of k y . Likewise, we assume 
that excitation and response of the system are monochro- 
matic [E,H ex exp(— iuit)]. The linearity of Maxwell's 
equations makes it possible to generalize the developed 
formalism for the arbitrary pulse excitation by repre- 
senting the incident radiation by linear combination of 
monochromatic waves. 

Note that in the process of calculating the waveguide 
modes, it may be necessary to determine the proper sign 
of the k x (or k z ) component of the wavevector in a partic- 
ular layer. If the component of the wavevector has com- 
plex value, this sign is determined from the requirement 



for the mode to be finite in its domain. If the wavevec- 
tor component is real, its sign should be determined to 
impose the propagation of energy in the positive x (or z) 
direction^. 

The set of waveguide modes defined above allows the 
introduction of the scalar product: 

/oo 
(Ei x Ml) ■ zdx (7) 
-oo 

where dagger (f ) corresponds to the adjoined field - field 
of the mode propagating in the reversed z directio n 13 ' 17 . 
It can be show n 13 ' 14 ' 17 that in a given multilayer (i) all 
TM-polarized waves are orthogonal to all TE-polarized 
waves, (ii) the guided modes are orthogonal to each other, 
and (iii) the top and bottom modes may have some cou- 
pling, depending on t\ and ejv: if ei = ejv, the top and 
bottom modes corresponding to the same value of k x are 
coupled to each other and are orthogonal to all other 
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modes; if one of the two materials is lossy (as it is usually 
the case with plasmonic structures), the top and bottom 
modes are, as a rule, orthogonal to each other. 

Note that similar to what has been suggested in RefJ^ 
the scalar product can be calculated analytically, signifi- 
cantly speeding up the calculation. 

III. MODE COUPLING ACROSS MULTILAYER 
STACKS 

A. General formalism 

We now turn to the main point of this work - discus- 
sion of coupling of the modes at the boundary between 
the two multilayer structures. For simplicity, we present 
results for the case when the interface is located at z = 
[FigQJc)]. Generalization of the technique for other lo- 
cations of the interface is straightforward. 

We are solving the classical scattering problem: find- 
ing the fields scattered by the interface provided that the 
incident fields are known. The incident fields are rep- 
resented by the modes propagating in the +z direction 
on the left-hand-side of the interface (z < 0) and by 
the modes travelling in the —z direction on the right- 
hand-side of the interface (z > 0); the scattered fields 
are represented by the modes travelling in the — z di- 
rection on the left-hand-side of the interface and by the 
modes travelling in the +z direction on the right-hand- 
side of the interface. The modal representation [Eqs.©] 
reduces the scattering problem to an arithmetic task of 
finding the coefficients and as a function of A\ 
and A^, which can be solved by imposing the following 
set of boundary conditions: 

El* = E Rm ; H Ly = H Ry , (8) 
El v = E Ry ; H Lx = H Rx . 

As can be explicitly verified, the remaining boundary 



conditions follow from Eqs.®. 

In the case of normal incidence (k y — 0), TM- and TE- 
modes do not couple to each other. Correspondingly, in 
this case the first two boundary conditions in Eqs.© 
describe the reflection, transmission, and scattering of 
TM-polarized waves, while the remaining two conditions 
describe the optical properties of TE-polarizcd modes. 

In order to solve the scattering problem, Eqs.© need 
to be converted into the set of coupled linear equations 
for the amplitudes of the scattered modes. To achieve 
this goal, we substitute the modal expansion [Eqs.ljSJ)] in 
Eqs.([8]), and subsequently multiply the resulting expres- 
sions by the adjoined fields of left- and right-hand-side 
modes, as illustrated below. 



B. Numerical implementation of the algorithm 

In numerical simulations, it is necessary to replace the 
continuous integration over k x with finite sums. Thus, 
Eqs.© becomes 

E = ^[A (m)+ E (m)+ + J 4 (m) "E (ra) > (ro) (9) 

II I 

H = ^[A (ro)+ H (ro)+ + A (m) ~H (m) "]t/; (m) 

m 

where E^" 1 ^ = E(fcj,, ±fci m ' > ), and similar for H; the sum- 
mation in Eqs.([9} goes over all modes (guided, top, and 
bottom), and the weight factors w are equal to 1 for 
the guided modes, and are determined by the integration 
method used for top and bottom modes^. Note that the 
number of modes on the left-hand-side of the interface 
does not necessarily equal to the number of modes on 
the right-hand-side of the interface. 
Eqs.© now become: 



JVx, N R 



\ W R ' 



m— 1 m— 1 



N L N R 



*R "it, \ W R 



ni—1 rn—1 



£ [4 m) 4'f + 4 mr < r ]4 m) = E t4 m) e£> + AfEjv (io) 



m— 1 m— 1 



JVx, N R 



To solve for N R + unknown amplitudes, we multiply first two equations in Eqs. (|10p by the fields of TM-polarized 
modes and integrate the resulting products over x; similarly, we multiply the two latter equations by the fields of 
TE-polarized modes and perform the integration. Assuming that the index m first spans the TE-polarized and then 
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TM-polarized waves, the procedure results in the following two sets of matrix equations: 
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M 4 



Si 



— ran . . 4 
+ A( m ) 
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(11) 
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£fi A R 



-£R 



ZR 



(m)~ 



(m)" 



(12) 



where the summation over repeated index m is assumed and matrix elements are given by: 



<R{L\R} 



9) 
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£.{L\R} 



>£,{L\R} 
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J {L|fl} 
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J {L|fl} 
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dx, n < iV^ E 



dx, n > N 
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{L\R}y 

o M {L|fl} x 
o M {L\R} y 



Rx 
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(n)" 

L y 

(„)- 



TE ' 



dec, n < iV^ E 
dx, n > 7V^ E 



dx, n < NT E 



(13) 



dx, n > N 
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TE ' 
R 



dx, n < TV™ 
dx, n > TV™ 



Note that the modes of z < region are not necessarily orthogonal to the modes in the z > region. Thus, the 



matrices ^^ L ,Sj 



across the interface. 



1 )KL ' ^£ft> 



S.R 



may have substantial non-diagonal components describing cross-talk of the modes 



In fact, the above matrices are square and invertible only when Nl = Nr, in which case one of Eqs. (|llll2p can 
provide the information required to solve the scattering problem. However, even in this case inversion procedure 
may lead to significant numerical problems, and is undesirable. When Nl ^ Nr these matrices are rectangular and 
thus, even theoretically, cannot be inverted. To overcome this difficulty, we reduce Eqs. (|llll2[) to the following set of 
equations that represent the generalization of transfer-matrix formalism for coupled waveguide structures: 



A R =^nAl+^A+ 

A+ =m 21 Ai+d\ 22 Ai 



A+ = £ 21 A R +£ 22 A+ 



(14) 
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(15) 



(16) 



Finally, we combine Eqs.(|14| and arrive to the generalization of the scattering-matrix technique that solves the 
problem of interlayer coupling: 



Al = (/-£i29t2ij £11 A- + (I-£ 12 m 21 ) £12^22 A+ 

A+ = (i -m^^Y 1 ^^! a r + (i-vhiZ^Y'm^ a+ 



(17) 



with / being identity matrices. Eqs. (|17p represent the 
main result of this work. 

Before illustrating the accuracy of the developed ap- 
proach, we would like to underline its main advantages 
and disadvantages. The approach is ideal for calcula- 
tion of light propagation in extended structures with rel- 
atively small number of multilayer segments; the increase 
of the number of segments results in additional memory 
use for each given segment and correspondingly, mini- 
mizes the advantages of wave-matching approach over 
finite-difference and finite-element schemes. 

The developed technique provides an effective mecha- 
nism for coupling the multilayers with high index con- 
trast (high optical mode density difference) by imple- 
menting multilayer-dependent number of modes. How- 
ever, our calculations show that careful design of spec- 
trum of the open- waveguide modes is necessary when the 
coupling to and from the highly-confined modes is calcu- 
lated. 

For guided modes, the approach allows straightfor- 
ward calculation of inter-mode cross-talk by calculating 
the mode-specific pointing-flux and multiplying it by the 
amplitude of the given mode squared; similarly, the ap- 
proach allow easy calculation of emission directionality, 
naturally separating the fields produced by guided modes 
from the fields of open-waveguide modes, and the fields of 
the incident waves from the fields of the scattered waves. 



IV. INTER-GUIDE COUPLING IN PLASMONIC 
AND METAMATERIAL SYSTEMS 

We now illustrate the accuracy of the presented ap- 
proach on several examples of plasmonic and metamate- 
rial systems. 

A. Light emission and scattering by a waveguide 

We first consider light coupling to and from the waveg- 
uide. As an example, we use a 600-nm-thick Si waveguide 
surrounded by air, and calculate the coupling between 
this system and a homogeneous dielectric at Ao = 1.5/zm; 
e Si = 12.1222. 

To analyze the accuracy of our technique we assume 
that the system is excited by TM 2 mode with ampli- 
tude of 1, and study the percentage of the reflected light 
into TMq,TMi and TM 2 modes as a function of dielec- 
tric permittivity of the homogeneous dielectric [see in- 
set in FigJ^a)]. For comparison, we have calculated the 
same parameters with commercial finite-element-method 
(FEM) software 2 ^. The perfect agreement between the 
results of our technique and FEM simulations is shown in 
Figj^l As expected even for this relatively simple system, 
the scattering-matrix approach uses orders-of-magnitude 
less memory than FEM model. More importantly, scat- 
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FIG. 2: (color online) (a) light reflection in the planar air- 
Si-air waveguide coupled to a homogeneous dielectric; Ao = 
1.5/^m, waveguide thickness d — O.Gfim; the geometry and 
profiles of the waveguide modes supported by the system are 
shown in the inset; excitation by TM2 mode is assumed; the 
graph shows the comparison between the presented here tech- 
nique (lines) and FEM simulations (dots); R2 = since the 
symmetry of TM\ mode is different from that of TMq and 
TM2 modes; (b) shows the reflection in plasmonic gap (Au- 
Air-Au) waveguide; Ao = d = 0.6/im; thickness of lines illus- 
trates the convergence of the computations; panels (c) and 
(d) illustrate the field distributions in the system (a) with 
e D = 12.12 (c) and eff = 3.6 + 0.05i, e y J = -12.2 + 1.36i (d) 



tering matrix formalism is easily extendable to the case 
of multiple guides or multiple interfaces. 

The field matching obtained with our technique for 
e D = 12.12 are illustrated in Figg] (c). Panel (d) of 
the same figure illustrates the matching obtained in cou- 
pling between an air-Si-air guide and an anisotropic hy- 
perbolic metamaterial with eff = 3.6 + 0.05z,e^ = 
-12.2 + LSQiS^^LSS^mi^^Si^. 



B. Light scattering in plasmonic systems 

To further analyze the accuracy and limitations of 
the developed field-matching technique we have calcu- 
lated the scattering from the plasmonic analog of the 
Si-guide presented above: a 600-nm-thick plasmonic gap 
waveguide operating at Ao = 0.6/im. We assume that 
cladding of this guide is composed from two gold strips 
with £Au = —8.94+ 1.32*22. The agreement between the 
mode matching and FEM calculations of modal reflectiv- 
ity in plasmonic gap guides is shown is Fig^b). 

We have further simulated the propagation of plas- 
monic modes across surface structure defects. In par- 
ticular, we have used our approach to simulate the SPP 
propagation in "plasmonic step" geometry (FigJ3|). As 
expected when the incoming SPP is travelling on the up- 



FIG. 3: (color online) Scattering of the SPP propagating 
across Au-Air step; panels (a) and (b) show transmission, 
reflection, and scattering (S — 1 — R — T) of an SPP that 
is incident on the step; geometry of the structures is shown 
in insets; lines correspond to the formalism developed in this 
work (line thickness represents data variation due to changes 
in spectra of bulk modes), dots represent results from FEM 
simulations; panels (c) and (d) illustrate the field distributions 
obtained from scattering matrix (c) and FEM (d) simulations 
for d — 0.7/im 

per side of the plasmonic step, most of the incident energy 
is converted into free-space modes. In contrast, when the 
incident SPP is travelling on the lower side of the step, 
the majority of energy is converted into the reflected SPP 
wave. 

Our numerical simulations demonstrate that at near- 
IR frequencies (when \e m \ ^> ed), the scattered field can 
be successfully decomposed into "top" open waveguide 
modes, as suggested in Ref. 1 ^.. However, in the proxim- 
ity of SPP resonance (e m ~ — e^), the inclusion of the 
"bottom" modes is necessary to adequately describe the 
optical properties of the system. 

Note that the results of wave-matching simulations 
are almost identical to those obtained with FEM. Once 
again, we underline that wave-matching technique allows 
for calculation of light propagation in much larger sys- 
tems than FEM system does. 

To assess the convergence of the our method, we have 
performed a set of simulations for each SPP structure de- 
scribed above, varying the configuration of spectra of top 
and bottom modes (for simplicity, equidistant k x spectra 
between 500 and 3000 modes were used). As expected, 
our simulations showed that it is necessary to design the 
spectrum of bulk modes to adequately resolve the SPP 
propagating at the z = interface. Interestingly, the 
"averaged" parameters (such as amplitudes of reflected 
guided modes) are much more sensitive to spectrum vari- 
ations than the matching of the boundary conditions at 
z = interface, which is often considered to be an indi- 
cation of accuracy of a numerical method. The typical 
inter-set variation of reflection, transmission, and scat- 
tering are illustrated in Figs l2l3l The latter figure also 



8 




0.04] 

< 

'< 

-0.04 



(c) 


Ah/ 












AT 








r 






T 


1 



0.02r 



(d) 



1000 2000 
mode number 



< 

-0.02- 



Ar 



TE.top "L 

TM.top 

TE.bottom 

TM, bottom 

1000 ' 2000 
mode number 



scattering or modal cross-talk. In this section we assume 
that the layers on both sides of the interface are aligned 
with each other (%L,j — x R,j)- 

In order to realize the efficient control over pulse prop- 
agation in the plane, the two layered structures must (i) 
have the same number of guided modes, and (ii) provide 
the ability to independently control the index of the mode 
[crucial for steering the light] and modal profile [crucial 
for optimizing the overlap integrals involved in £ and 
matrices] . 

As shown in Rcf. 1 - on the example of surface waves, 
these conditions can be satisfied when 



~2 xx 

n e T . 



(18) 



FIG. 4: (color online) (a,c) an interface between air-Si-air 
waveguide and isotropic air-(e = 6.06)-air guide leads to sub- 
stantial modal cross-talk, polarization mixing, and out-of- 
plane scattering; while the interface between air-Si-air waveg- 
uide and its anisotropic truly planar optics analog allows for 
ideal mode matching with light steering capabilities (b,d); 
guided modes in (c,d) correspond to mode number > 2000; 
the system is excited by a TM2 guided mode propagating at 
the angle 30° to the z — interface; the amplitudes of modes 
in panels (c,d) are normalized to the amplitude of the incident 
mode 



shows the agreement between the field distribution ob- 
tained with our generalized scattering matrix formalism 
and with FEM. 



V. TRULY PLANAR OPTICS 

We now turn to the analysis of inter- mode coupling and 
out-of-mode scattering in planar optics. Planar optics in 
general (and SPP optics in particular) are fundamentally 
different from their free-space counterparts. Thus, when 
a plane wave is incident on the planar interface, the scat- 
tered field can be decomposed into one reflected plane 
wave, and one transmitted plane wave. In contrast to 
this behavior, when one guided mode is incident on the 
planar interface between two waveguides, it generates a 
continuum of open-waveguide modes in addition to the 
(sets of) reflected and transmitted guided waves. 

As seen from Figs l4l5l and from RefJ^, the typical in- 
terface between waveguide systems leads to scattering 
of about 20% of incident radiation. Every attempt to 
change the effective index of the mode in planar struc- 
tures is necessarily accompanied by modal cross-talk or 
by out-of-plane scattering of radiation. 

It is possible, however, to utilize anisotropic metama- 
terials to completely eliminate the cross-talk, and to map 
the familiar laws of 3D optics to optics of planar guides. 
The main idea of planar optics lies in the ability to guide 
light along the planar optical circuit with no out-of-plane 



with n being the constant number that does not depend 
on layer number j. 

As can be explicitly verified, when Eq. (fl~g|) is satisfied, 
all € and $3 matrices become diagonal. Thus, the inter- 
mode coupling is absent across the interface. The in- 
terface remains completely transparent to TE-polarized 
waves, while reflection and refraction of TM-polarized 
modes are controlled by the ratio of out-of-plane permit- 
tivities. The direction and amplitudes of the reflected 
and refracted modes are related to the direction and am- 
plitude of the incident modes via the following Snell's 
law: 



sin [Vi — sin 



and Fresnel coefficients: 
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(20) 



The above equations represent generalization of the for- 
malism of RefJ£ to multilayered guides. 

The concept of truly 2D optics is illustrated in FigfJ] 
on example of air-Si-air system coupled to metamateri- 
als waveguide. As expected, reflection of a single mode 
in conventional planar system is accompanied not only 
by significant radiation scattering and modal cross-talk, 
but also by cross-polarization coupling. In contrast, for 
metamaterial guides the single incident mode excites one 
reflected wave, and one transmitted wave. 



A. Low-scattering plasmonics 

An interesting extension of the truly planar optics is 
possible in plasmonic systems. From the fabricational 
standpoint, it is highly desirable that the plasmonic cir- 
cuit is fabricated on top of common metallic substrate. 

Fabrication of extremely low-scattering plasmonic cir- 
cuit is possible when using anisotropic dielectrics de- 
posited on noble metals in the limit of visible and 
near-IR frequencies, where the permittivity of metal is 



9 




approximately given by 



2 
z, um 



FIG. 5: (color online) Anisotropic coatings (b,d) can signifi- 
cantly reduce (and almost eliminate) the scattering losses in 
plasmonic circuits; the figure shows reflection, transmission, 
and scattering in conventional plasmonic circuit (a); field dis- 
tribution for e_D = 4 (c) clearly shows that only a fraction 
of the energy of incident SPP is transferred into SPP on the 
right-hand-side of the interface; on the other hand, the inter- 
face between anisotropic system with e yz = 1, e M = 4 allows 
for substantial modulation of SPP index, and results in almost 
perfect SPP-SPP coupling across z = interface 



much larger than the permittivity of dielectric (\e m \ 3> 
e d i e d )• 

Consider the situation where the in-plane (yz) compo- 
nent of permittivity of dielectric is kept constant across 
plasmonic system, and only out-of-plane (xx) component 
of permittivity of the superstrate is varied. In such a 
system the propagation constant of the plasmonic mode, 



!J~- 



(21) 



can be effectively controlled by changing the parameter 



In 

cay 



the 
of 



same 
the 



limit, 
mode 



the exponential de- 
into the dielectric 



does 



E,H cx exp(-Kdx), n d ~ w/cy e y d z /\e m \ 

not depend on its propagation constant. The 

only source of out-of-plane scattering in such a 

structure is related to a weak dependence of the 

field profile in metal [E,H cx exp(K m x)] on e^ x : 

n m ~ w/<V|em|(l - e^/kml)- 

As seen from FigJSJ anisotropy of dielectric super- 
strate, achievable, for example, with electrooptic effect—, 
provides orders-of-magnitude suppression of out-of-plane 
scattering in plasmonic systems with respect to isotropic 
counterparts. 

VI. CONCLUSIONS 



To conclude, we have developed a reliable numerical 
technique for calculation of light propagation in planar 
guides and in arrays of planar guides. We have illus- 
trated the developed formalism on examples of photonic, 
metamaterial, and plasmonic guides, and presented an 
approach to utilize anisotropic metamaterials for mini- 
mization and elimination of modal cross-talk in planar 
optical circuits. 

This work was sponsored by the ONR (Grant No. 
N00014-07-1-0457), NSF (Grant No. ECCS-0724763), 
and AFOSR (Grant No. FA9550-09-1-0029). 
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